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• A fireball model with time evolution based on transport calculations is used to 

, examine the dilepton emission rate of an ultra-relativistic heavy-ion collision. A tran- 

' sition from hadronic matter to a quark-gluon plasma at a critical temperature Tq 

OO \ between 130-170 MeV is assumed. We also consider a possible mixed phase scenario. 

■ We include thermal corrections to the hadronic spectra below Tc and use perturbation 

theory above Tq. The sensitivity of the spectra with respect to the freeze-out tem- 
perature, the initial fireball temperature and the critical temperature is investigated. 



Abstract 



o 
o 

p ,' 1 Introduction 

(D ■ 

' There are convincing arguments that the theory of the strong interaction, QCD, exhibits 

a phase transition at sufficiently high temperatures from a confined hadronic phase to the 
^ . quark-gluon plasma (QGP) phase Lattice QCD results set the critical tempera- 

^ ; ture Tc to about 150-200 MeV [|] |][|]. It is hoped that it is possible to create the QGP 

in ultrarelativistic heavy- ion collisions at CERN and BNL, and there are preliminary in- 
dications that the QGP might have already been encountered 

Dileptons {e^e" or fi^fi^ pairs) are considered to be an optimal probe for the early stages 
of the collision because they leave the hot region without thermalization. Hadrons which 
reach the detector can only tell us about the later stage of the collision, the freeze-out 
zone. When the fireball created in a heavy-ion collision is in the QGP or partonic phase, 
dileptons are mainly produced by thermal quark-antiquark annihilation. As the system ex- 
pands, it cools off and undergoes the transition into the hadronic phase. There, dileptons 
come mainly from pion and kaon annihilation processes, which are dynamically enhanced 
through the formation of light vector meson resonances (p, to and (p with masses below 1.1 
GeV). The dilepton emission rate is proportional to the imaginary part of the hadronic 
current-current correlation function, or the electromagnetic spectral function. The in- 
variant mass distribution of the lepton pairs reflects the mass distributions of the vector 
mesons at the instant of decay. This offers the possibility to study the influence of finite 
temperature and baryonic density on the meson spectral distributions. 
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In this paper, we calculate the purely thermal dilepton emission rate for both phases and 
compare it with data from the CERES /NA45 experiment. Much work has been done 
to interpret the observed dilepton enhancement in the low-mass region BP] in terms of 



in-medium modifications of hadron spectra |24]. We investigate the thermal modi- 



fications that the hadrons may exhibit in the limit of zero baryon density. We demonstrate 
that these alone can also explain the data very easily if the initial fireball temperature 
is sufficiently high to support a QGP during at least part of the expansion. Finally, we 
study the dependence of the shape of the dilepton emission spectra on the parameters of 
our simple fireball model and discuss a recently proposed mixed phase structure. 

2 Dileptons from a fireball 

The dilepton emission rate from a hot domain populated by particles in thermal equilib- 
rium at temperature T is proportional to the imaginary part of the spin-averaged, time-like 
photon self-energy, with these particles as intermediate states. The thermally excited par- 
ticles annihilate to yield a time-like virtual photon with four-momentum q which decays 
subsequently into a lepton-antilepton pair. The differential rate is given by 

where a = e^/An, (3 = 1/T, and we have neglected the lepton masses. We defined 
ti{q) = — Here Yl^^ denotes the trace over the thermal photon self-energy which is 
equivalent to the thermal current-current correlation function 

n^,(g,T) = ijd^x e'''^{Tj^{x)M0))i3, (2) 

where is the electromagnetic current. The result eq.(|l|) is valid to order a in the 
electromagnetic interaction and to all orders in the strong interaction. 

To compare with experimental data, we set up a model for the space-time evolution of a 
heavy-ion collision, assuming approximate thermal equilibrium to be a useful concept pi 



Some recent discussions suggest |13][14] that equilibration times at the conditions of heavy 
ion collisions at CERN may indeed be very short (less than a fermi), small compared to 
expansion times of order 10 fm/c, although this is still a matter of debate. The easiest 
approach for our purposes is to use a simplified fireball model which parametrizes the time 
dependence of temperature and volume in accord with microscopic transport calculations 



[16|[17|. This assumes that the fireball can be characterized by a homogeneous temperature 



at all times and that it cools off adiabatically with the parametrization 

T{t) = {T - T°°)e-*/^i + r°°. (3) 

This ansatz introduces the initial temperature T*, a time constant ri and an asymptotic 
temperature r°°. Furthermore, we assume an isotropic expansion of the fireball in the 
centre-of-mass frame of the collision such that the time evolution of the volume can be 
described by 

ni) = ^, p(t) = A-*/-, (4) 

where Nb is the number of baryons which participate in the reaction, T2 is another time 
constant and is the initial baryon density of the hot spot. More sophisticated volume 
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parametrizations can be used, but they introduce more parameters and the differences in 
the resulting integrated rates are marginal. 

Finally, to obtain the measured dilepton rates, we integrate eq.(||) over the space-time 
history of the collision to compare them with the CERES/NA45 data. The CERES 
experiment is a fixed-target experiment, i.e. the dilepton rates are measured in the lab 
frame whereas our volume parametrization (Q) is valid for a fireball at rest, in the cm. 
frame of the collision. In the lab frame, the CERES detector covers a limited rapidity 
interval r/ = 2.1 — 2.65. The rapidity needed in order to boost to the cm. frame is therefore 
2.375. The CERES collaboration chose to display its data in the format d?N/dMdr], where 
rj is the (pseudo-)rapidity of the virtual photon. The longitudinal momentum pL is related 



to rj hy PL = niT sinhr] with the 'transverse mass' = \J M"^ + Pr- To be able to 
compare our calculated rates to the data, we integrate these rates over the transverse 
momentum px only, given that 

ct^p = MpT dM d7] dpT de. 

This is equivalent to integrating over all three-momenta as long as we set 77 = in the 
cm. frame. 

The formula for the space-time- and pr-integrated dilepton rates is now 



d^N f r 



dN{T{t),M,r,,PT 



d'^xd'^p 



Acc{M,r],pT), (5) 



where t/ is the freeze-out time of the collision, and the matrix Acc{M,rj,pT) accounts for 
the experimental acceptance cuts specific to the detector. At the CERES experiment, each 
electron/positron track is required to have a transverse momentum pT > 0.2 GeV, to fall 
into the rapidity interval 2.1 < rj < 2.65 in the lab frame and to have a pair opening angle 
@ee > 35 mrad. Finally, for comparison with the CERES data, eq.(^) has to be divided 
by dN/dr], the number of particles per unit rapidity. As the data have been taken close 
to midrapidity, i.e. around the central plateau of the approximate Gaussian distribution 
of dN/dr] in the cm. frame as mentioned above, its value is assumed to be constant over 
the r] range restricted by the detector cuts. 



3 Spectrum above critical temperature 

Lattice QCD calculations suggest [|| [^J that sufficiently far above the critical temperature 
Tq for the deconfinement transition, quarks and gluons can be treated approximately as 
an ideal gas; they interact only weakly because of the behaviour of the running coupling 
strength as at high T. In this temperature region we can evaluate the hadronic part of 
the photon self-energy using perturbative QCD. To lowest order in a and to zeroth order 
in as, n^jjy corresponds simply to the thermal quark-antiquark loop to which we restrict 
ourselves in the present paper. Higher order thermal QCD corrections beyond the leading 
qq loop will be investigated in forthcoming work. 

At our typical temperatures T of several hundred MeV, we need to take into account the 
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three flavours u, d and s. Heavier quarks such as the c and b are strongly suppressed 
by Boltzmann factors because of their large mass and can be considered as 'frozen'. The 
total contribution to the dilepton emission rate is obtained by summing the individual 
qq contributions multiplied by the squares of the electric charges e/ of their respective 
flavour currents. In addition, each quark exists in three colours which leads to an overall 
multiplicative factor of Nc = 3. Furthermore, the u and d quarks are considered massless 
and the heavier ss-pair can only be produced if q'^ is larger than 4mg. So the total 
imaginary part of the self-energy which will be used in eq.(|T|), evaluated with standard 
flnite-temperature techniques |1S], is 



lmU{q,T) 



where 



and 




fniE) 




(6) 
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In the limit T — > 0, eq.(^) reduces to the well-known vacuum result. At flnite T, the 
quantity (^) is always smaller than Imn(g,T = 0) as a consequence of Pauli blocking. 



4 Spectrum below critical temperature 

At T < 150 MeV, conflnement sets in and the effective degrees of freedom of the QCD 
Hamiltonian are now colourless hadrons. The photon can directly couple to = 1~ states 
(the lowest "dipole" excitations of the QCD vacuum). In the energy region of interest, 
these particles are the p, oj and (p mesons and multi-pion states carrying the same quantum 
numbers. Our next step is thus to connect the electromagnetic current-current correlation 
function with the currents generated by these mesons. We use an effective Lagrangian 
which approximates the SU{3) flavour sector of QCD at low energies. The appropriate 
model for our purposes is the improved Vector Meson Dominance model combined with 



chiral dynamics of pions and kaons as described in [19| 



Within this model, the following relation between the imaginary part of the irreducible 
photon self-energy Imll and the vector meson self-energies Ilv{q) in vacuum can be derived: 

Imn(,) = (T) 

V 3v 

9 \ ^2 



— rriy + iimliv{p) 
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where my are the (renormahzed) vector meson masses Q gy is the couphng and g 
the couphng, where ^ stands for one of the pseudoscalar Goldstone bosons 7r^,7r'^ 
and K^,K^. Eq.(j^ is vahd for a virtual photon with vanishing three- momentum, q = 
(g*^, 0,0,0). Taking the hmit |g| = should be reasonable for our purposes in view of 
the fact that the cm. rapidity interval accessible at CERES restricts |g| on average to 
a fraction of the vector meson mass my- For finite three-momenta there would be two 
scalar functions tl^ and IIt, because the existence of a preferred frame of reference (the 
heat bath) breaks Lorentz invariance. In the following, we calculate the hadronic spectra 
liy at finite temperature. 



4.1 The p meson spectrum 

The main decay channel of the p-meson is /? ^ vr+vr" with a vacuum width V ~ 150 MeV. 
In the VMD model, hadronic current conservation leads to two possibilities to couple the 
p to pions: There is a pion tadpole diagram and a two pion loop diagram. Both terms 
contribute to lip and correspondingly modify the p meson properties. Their vacuum values 
are well-known [|l9|, thus in the following we will calculate only the changes induced by 
finite temperature. 



Using the real-time formalism of thermal field theory [18|, we find for the imaginary part 
of the two pion loop diagram 



Imn(,°, q = 0) = (l - ' (l + VB{q'm) 0{q' - 2m^). (9) 

This is exactly the T = result multiplied by the characteristic thermal factors: the imag- 
inary part is larger than at zero temperature as a consequence of the Bose-enhancement of 
the thermal pions. Using the relation Imll(mp) = nipT, we see that the decay width Tp-^TrTr 
of the p meson is enhanced because there are already thermally excited pions present in 
the heat bath. At temperatures of about 100 MeV, this effect causes an increase of the p 
width by about 20% . Additional sources of thermal broadening at high temperatures are 
reactions of the p meson with pions in the heat bath, such as pvr — > Svr. We checked that 
this is only a small correction to the width. Other decay modes of the p like p — 7r"'"7r~7, 
not taken into account, may also start to play a role at high T, however, their branching 
ratios are very small. 

The calculation of the real part of the thermal p meson self-energy is technically more 
involved, but the resulting temperature dependent mass shift 5mp{T) is small. Eletsky 



and loffe [ 20 1 have pointed out that the leading term of 6m p starts at order when 
p — ai mixing is properly taken into account. The actual evaluation of 5mp gives an (up- 
ward) mass shift of less than 3% of mp even at temperatures as high as 150 MeV. We can 
therefore safely neglect this effect in our calculation. 

Note that the general form of the thermal p spectrum does not change dramatically, at 
least up to T < 150 MeV. Substantial broadening of the p meson spectrum is expected 
in the presence of large baryon density |l^|24]. However, under the conditions of Pb-Au 



collisions at CERN, the fireball rapidly expands and the baryon density at hadronization 



^In the present context, mv is understood to include the shift from bare to physical mass induced by 
ReHv PI. 
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drops to a fraction of normal nuclear matter density. We can thus ignore such density 
effects. 



4.2 The uj meson spectrum 

The main decay mode of the a;(783) vector meson with a branching ratio of 89% is w — > 
TT+vr'^Tr" and T^^^^■,^ = 7.5 MeV. It is well-known that a point-like coupling a; ^ Svr 
does not describe the oj decay properly. It is necessary to include an intermediate step, 
the Gell-Mann, Sharp, Wagner (GSW) process cj ^ pvr ^ 7r+7r~7r''. The full decay 
process is described by an interference between the amplitudes of these two processes. In 
the following, we assume that the intermediate virtual rho meson propagation is merely a 
'vertex' correction and that the effects of temperature on the p meson are small, so we can 
neglect them to leading order. This is justified by noting that the thermal production of a 
p meson is heavily suppressed by e~'^pl'^ ^ 1 in the temperature range we are interested 
in (up to about 150 MeV). 



The explicit calculation ||2^ for an uj meson at rest, = (w, ^= 0), leads us to the result 

Imn(a;) = ^^(^i+SiSs), (10) 

Si = I dE+dE^ {q^q^ - {q+ ■ q^f) \M{io; q+, qo, q-)\^ ■ n{E+, E^,u;), 
Ai 

B2 = - J dE+dE^ [q^q^ - {q+ ■ q^f) |A1(c^; g+, -go, Q-)P • n{E+,E^,uj), 



where 



\AI \ , o 3VVP 

M{uj; q+, qo, q-) = —pr- + 2 — - — gp. 



/I fn a=+fi,^ (^^ ~ -mj- Up{q^ - qa) ' 

Here, Hp is the zero-temperature rho meson self-energy. Its explicit expression and a 
general discussion of this matrix element can be found in [^9|, including the coupling 
constants h and gwp ifn is the pion decay constant). E± and q± are the energies and 
momenta of the charged pions in the final state. To keep our notation short, we have 
introduced the function 



n{E+,E_,u;) = ((1 + /B(i?+))(1 + /B(i?-))(1 + /b(^ - - i^-)) 
-fB{E+) fB{E-) feiuj -E+-E^ 



In eg. (p^) , Bi describes the process of a Bose-enhanced three-particle decay of the oj me- 
son, n{E^,E-,uj) accounts for the characteristic Bose factors. The kinematically allowed 
integration region in the E^E- plane, Ai, is the same as in the T = case for the decay 
process. It is limited to rriT^ < E± < uj — 2171-^ and 2/71,^ < -E+ + -E- < uj — m-,^ with the 
constraint q^q3 — {q+ ■ q-Y > 0- 

To interpret the B2 term in eq.(lO), we note that 

-n{E+,E_,uj) = [{1 + fB{E+mi + fB{E^)) fB{E+ + E, - uj) - 

-fB{E+) fB{E-) (1 + fB{E+ + E^- u:))) , (11) 
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by using 1+ fB{E) + fs^—E) = 0. Thus B2 corresponds to the scattermg of the uj meson off 
a thermally excited vr*^ into tt^ and vr" , where each particle has the characteristic thermal 
corrections. The integration region A2 is still constrained by q^q3 — {q+ • Q-)'^ > 0, but 
now in the region defined by + E- > a; + m-,^. Because the integration region A2 is 
now unbounded and the functions involved are only damped by Bose factors and their 
combinations, the lott — > tttt scattering contribution to the decay width increases quite 
strongly with temperature, as we will see below. The factor 3 multiplying B2 comes from 
summing over all possible charge combinations co + 7r~^ ^ tt^ + n'^ , uj + ^ + 7t~ and 
(J + vr*' — > vr^ + 7r+. 

The evaluation of the real part of the lv meson self-energy coupled to the Svr continuum 
is much more difficult than that of the imaginary part because it would involve a two- 
dimensional Cauchy principal value integral. We refrain from calculating the real part of 
the three-body decay, assuming that its temperature dependent mass shift is not much 
different from the (small) one of the p meson. 

4.3 The (f) meson spectrum 

Because of the large mass of the (j) and its main two-body decays into heavy kaons, the 
impact of temperature on its spectrum is marginal. The thermal production of kaons 
at temperatures even as high as 150 MeV is still suppressed by an order of magnitude 
compared to thermal pion production. We find that the sharp peak structure of the (p 
meson does not change as its total width increases by about 20% to just 5.2 MeV at 
the highest temperatures. Thus it is only moderately affected by temperature effects 
and one should be able to observe its sharp resonance structure in future dilepton rate 
measurements. 

5 Integrated rates 

Whereas the p and (j) mesons experience only moderate thermal changes of their spectral 
functions, the oj spectrum changes its shape considerably, mainly due to ujtt tttt reactions 
in the heat bath. We now use eq.(0) to relate our temperature-dependent vector meson 
self-energies to the imaginary part of the electromagnetic current-current correlator which 
enters the dilepton rate formula (|l|) as long as we stay below Tc- Above this critical 
temperature, we use expression (^) which yields the dilepton emission rate from a thermal 
quark-antiquark system. It is instructive to have a look at the spin-averaged spectral 
distributions R{q,T) = l2TTlmIl{q,T) /q'^ at different temperatures T (presented in fig.^ 
for q^^ = {u},q = 0)). Obviously, the spectral functions on either side of the quark-hadron 
transition are quite distinct: the qq loop spectrum is flat with no significant structure 
whereas the hadronic spectra exhibit a clear resonance structure. How these differences 
manifest in the dilepton emission spectra is examined in this section where we convolute 
the dilepton rate with the space-time evolution of the fireball system to obtain the total 
rate (^) which can then be compared to the CERES/NA45 data. 

For a start, we use the parameters shown in table || for eqs.(|3|) and (^) which are adapted 
to reproduce the results of microscopic transport calculations jl6|||17[| for the space-time 
evolution of the S-Au and Pb-Au collision (we show the Pb-Au parameters as this is the 
case of primary interest). Let us now apply our "mixed" model scenario (hadronic degrees 
of freedom below Tc, quarks above Tc) with Tc = 150 MeV and initial temperatures 
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Figure 1: Spectral distribution R = 127rlmn/g^ for = {u},q = 0) for different tempera- 
tures above and below the quark-hadron transition. The background contribution to the 
hadronic spectra at large u arises from the 47r continuum. 



of 190 MeV and 210 MeV, respectively. These two cases are shown in fig-^. For the 
moment, we employ a direct transition between the two phases at Tc and comment on the 
implications of a possible mixed phase later. We find that we can reproduce the shape of 
the data quite well. We need a normalization factor k which can be explained as follows: 
for the hadronic phase where thermal tttt annihilation is dominant, there is an increased 
population of the pion phase space due to a finite pion chemical potential we did not include 



[24|. However, the striking feature is that when choosing a higher initial temperature, 
T* = 210 MeV, the normalization factor needed to reproduce the data becomes smaller, 
consistent with the picture that the system now stays for almost the whole expansion time 
in the partonic phase where there is no pionic chemical potential. We also checked that 
our calculations reproduce the observed transverse momentum distributions: the main 
part of the enhancement resides in the px < 500 MeV region. 

We can therefore easily account for the dilepton "excess" in the intermediate mass region. 
At high masses, the scenario with T* = 190 MeV still underestimates the data a bit. With 
a higher initial temperature T* = 210 MeV, the slope becomes flatter and moves up in the 
right direction. This can be understood as follows: using a spin-averaged spectral function 
as in fig. [l|, we can rewrite eq.([l|) as 

dN a2 R{q^,q) 



d^xd^q Utt^ ePi'' - 1 ' 
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From eq.(|^) it follows that the thermal spectrum R approaches its constant perturbative 
plateau for large momenta and high temperatures or for any momentum at low tempera- 
tures. From this we infer that in the QGP phase, the shape of the dilepton emission rate 
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Pb-Au 



rpi 


190 MeV 




105 MeV 


n 


10 fm 


p' 


2.55 po 


T2 


6 fm 


Nb 


260 


dN/dri 


220 


tf.o. 


10 fm 


Tf.o. 


136 MeV 


Pf.o. 


0.48 po 



Table 1: Typical set of parameters controlling the space-time evolution of the fireball, 
adapted for the Pb-Au collisions at CERN-SPS. Baryon densities are given in units of po 
= 0.17 fm-^ 



for high invariant mass at high T is basically determined by the inverse of the temperature, 

/3, as 

^ ~ const xe-*. (13) 

As the integrated rate (|5|) averages over all temperatures according to the fireball 
parametrization, we see in fig.§ a straight fall-off at high invariant masses, with a slope 
parameter — (1/T), averaged over the temperature evolution of the fireball. This is an 
important result: if we start with a different parametrization of the fireball, say, with a 
higher initial temperature, the slope changes, it becomes flatter because we average over 
1/T. Thus an accurate measurement of the dilepton rate at high invariant mass yields 
valuable information about the temperature evolution of the fireball, assuming that equi- 
libration is fast. 

In contrast to its appearance as a sharp resonance at T = 0, the uj meson is hardly visible 
as a little bump on top of the broad p meson spectrum. This was to be expected because 
the temperatures of the system in the hadronic phase before freeze-out are still quite high, 
between 150 and 125 MeV. 

The (p meson can still be identified as a sharp resonance. As the (p is hardly influenced by 
temperature, its height over background ratio may thus act as a chronometer for the criti- 
cal time tc once the fireball parameters are sufficiently well known: at T* = 210 MeV, the 
(j) meson still sticks out of the QGP 'background', but with a reduced height, confirming 
that the system stays now longer in the QGP phase than in the scenario with T* = 190 
MeV. 

To get a better understanding for the space-time evolution of the fireball, let us look at the 
time snapshot picture of the fireball expansion in fig.^. Shown are the integrated dilepton 
rates between the start of the expansion and the freeze-out time tf = 10 fm/c in steps 
of At = 0.5 fm/c. Although the QGP phase is very hot in the beginning (to = 0), its 
contribution to the dilepton rate is initially small because it occupies only a small volume. 
With the consecutive expansion, the rate per unit volume decreases as the temperature 
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Figure 2: A mixed scenario for the Pb-Au case with two different initial temperatures. 
The background at low energies consists mainly of vr^, r/ and uj Dalitz decays. The data 
are taken from Pl. 
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Figure 3: A time snapshot of the Pb-Au fireball expansion with T* = 190 MeV (no k 
factor). Shown are the integrated dilepton rates between the start of the expansion at 
t = and the freeze-out time, tf = 10 fm/c in steps of At =0.5 fm/c. The CERES 
acceptance is included. The different behaviour of the partonic and the hadronic phases 
below and above Tc = 150 MeV is clearly visible. 
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Figure 4: The resulting Pb-Au dilepton emission spectra for three different critical temper- 
atures, starting from the same initial temperature T* = 210 MeV. For the sake of lucidity, 
the low-mass hadronic background has been omitted. No k factor has been applied. 



goes down. However, the increase in volume outweighs this effect partly. The shape of 



the spectrum stays roughly the same, but the slope for high M is now, following eq.(13), 

—I3\J -|- integrated over the transverse momentum px- The change in the slope with 
increasing time, or equivalently, with decreasing temperature is small yet visible, and the 
magnitude of the slope is set by the initial temperature. 

At tc" = 7 fm/c, the phase transition occurs, the system goes into the hadronic phase and 
the p and (j) meson resonances start to grow out of the QGP 'background'. It is clearly 
seen that the shape of the dilepton spectrum at high masses M is only determined by the 
QGP phase and that the 4> meson peak height would be larger if the transition occurred 
earlier or the freeze-out time were later. 

With these observations in mind, we investigate the sensitivity of the spectra with respect 
to the parameters Tc and tf. In fig.0, the resulting Pb-Au dilepton emission spectra are 
shown for three different critical temperatures Tc =130, 150 and 170 MeV. The initial 
temperature is 210 MeV. 

The slopes for high invariant mass M are equal in all three cases of Tc, as expected. For 
the lowest transition temperature, Tq = 130 MeV, we obtain a pure QGP spectrum. The 
higher Tc is chosen, the more 'hadronic' the spectrum looks and the less dileptons are 
present in the intermediate mass range 0.3 — 0.6 GeV. At high Tc, the broad p meson res- 
onance sticks out clearly, so does the narrow (p meson. From its peak height it is possible 
to judge how long the system has stayed in the hadronic phase before freeze-out. Note 
that the uj meson is not visible even for high Tc, because of its increased thermal width. 
An amusing side note is that all three curves have four points in common, left and right 
to the resonances. We also checked that the shape of the dilepton rate is not sensitive to 
the freeze-out time. 
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Figure 5: Fireball temperature as a function of time for a direct transition and a possible 
mixed scenario. 



So far, we have assumed that the transition from the QGP to the hadronic phase is instan- 
taneous, in the sense that the quark phase at T > Tq is separated from the hadron phase 
at T < Tc, and does not affect the time evolution of temperature and volume. Exactly 
how the hadronization occurs is not known at this stage, but there are indications that 
it may actually extend over a finite time interval, during which a mixed phase prevails 



[25|[26|. Lattice data show a discontinuity in the energy density at Tc so this stage is 
accompanied by a release of latent heat. Such processes can imply that the temperature 
stays almost constant during the transition time which may be considerably long, of the 



order a few fm/c |2^|28]. In order to estimate the consequences of these effects on the 
dilepton spectra, we have modified the time evolution of the fireball temperature such that 
there is a plateau of constant temperature Tc, as shown in fig.^ The width of the plateau 
is denoted by At, and the initial and freeze-out temperature are kept at the values listed 
in table During the transition, we assume that both phases exist in a homogeneous 
mixture and their relative occupation of the volume changes with time. Correspondingly, 
the spectral function entering the dilepton rates in eq.([l2D is modified to 



Rit) = a{t)RQGp{Tc) + [1 - a{t)]Rhadr{Tc), (14) 

where a{tc — At/2) = 1, a{tc + At/2) = and between these limits a drops linearly. We 
found that even if At is set to a value as large as 5 fm/c, the difference in the dilepton 
rates is at most 40%. For At = 1 fm/c, the rates are virtually indistinguishable. Within 
the error bars of the data, we conclude that the existence of a mixed phase is neither 
confirmed nor excluded at the resolution currently available. 
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6 Conclusions 



Within our simple fireball model, we can indeed reproduce the observed CERES dilepton 
spectra within a pure thermal model, assuming that there exists a partonic phase above 
a certain temperature Tc- The transverse momentum dependence is correctly described, 
too. With a sufRcicntly high initial temperature = 210 McV, the shape of the spectrum 
is perfectly well described. At such high initial temperatures, the system stays in the quark 
phase for much of its expansion period, and the dilepton radiation originates primarily 
from quark-antiquark annihilation, resulting in a flat spectrum. Hadronization in the final 
stage then leads to the appearance of resonance structures. We note that the thermal life 
time of the uj meson may decrease such that the u) actually decays inside the fireball. The 
detailed shape of the spectra is sensitive to certain combinations of freeze-out, initial and 
critical temperature. 

An improved experimental mass resolution in the region M > 1 GeV may yield important 
information about the initial temperature of the partonic phase. The role of the (f) meson 
as both a thermometer and a chronometer for the fireball expansion has been emphasized. 
Finally we investigated a possible mixed phase scenario and found that it would not leave 
any distinct traces in the shape of the dilepton rates. 
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